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Abstract 

We construct a family of Fourier Integral Operators, defined for arbitrary large times, 
representing a global parametrix for the Schrodinger propagator when the potential is 
quadratic at infinity. This construction is based on the geometric approach to the corre- 
sponding Hamilton-Jacobi equation and thus sidesteps the problem of the caustics gen- 
erated by the classical flow. Moreover, a detailed study of the real phase function allows 
us to recover a WKB semiclassical approximation which necessarily involves the multi- 
valuedness of the graph of the Hamiltonian flow past the caustics. 
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1 Introduction and statement of the results 

Let us consider the initial value problem for the Schrodinger equation: 

( 

ihdtipit, x) = -—ATp{t, x) + V{x)ilj{t, x), 

2m (^^1) 

'0(0,2;) = ip{x). 

where the potential V £ C°^(M";R) is assumed quadratic at infinity. In this case it is well 

known that the operator H in defined by the maximal action of A + V(x) is 

2m 

self-adjoint. Hence the Cauchy problem (jl.ip considered in L^(M") admits the unique global 
solution iP{t,x) = e-^^*/^(/7(x), Vt G M, V99 G ^^(M"). 

Under the present conditions a parametrix of the propagator under the form of a semi- 
classical Fourier integral operator (WKB representation) has been constructed long ago by 
Chazarain [Ch| (for related results by the same technique see also [Fuj . |Kij : for recent re- 
lated work see |MYj , [Yalj , |Ya2j ) . The occurrence of caustics of the Hamilton- Jacobi equation 
makes this construction local in time; the solution at an arbitrary time T > requires multi- 
ple compositions of the local representations. A global parametrix for the propagator has been 
constructed through the method of complex valued phase functions (as in }KSj . |LSj ). with 
related complex transport coefficients. A particularly convenient choice of the complex phase 
function (the Herman-Kluk representation) has been isolated in the chemical physics litera- 
ture long ago (} H-K ]). Its validity has been recently proved in |SwRj and |Ro2j ). The relation 
between the above approaches and the underlying classical flow is however less direct than 
the standard WKB approximation in which the phase function solves the Hamilton-Jacobi 
equation. 

In this paper we study the problem through the geometric approach to the Hamilton-Jacobi 
equation (see e.g. |CZ] . [Sik86| ). In Theorem 1.1 a parametrix is obtained for the propagator 
U{t) := e*^*/^ valid for t £ [0,T], < T < 00, under the form of a family of semiclassical 
^r/o&a/ Fourier Integral Operators (FIO), which extend to continuous operators in L^(M"). The 
corresponding phase function is real and generates the graph of the flow of the Hamiltonian 

T-L = h V(x). This technique not only yields globality in time, but also helps to obtain a 

2m 

unifled view of Fujiwara's as well as Chazarain's approaches on one side, and of the Laptev- 
Sigal one on the other. In Theorem 1.2 we prove that a WKB construction is still valid, 
necessarily multivalued because of the caustics. 



We assume: 

V{x) = {Lx,x) + Vq{x), L£GL{n), det L ^ 0; (1.2) 
Vo G C°^(M"), \d^Voix)\ < Co. (1.3) 

Then the main result of the paper is: 

Theorem 1.1. Let i fTTl)) and i fTT^) be fulfilled. Let < T < oo and ip £ S{W). Then: 

il>{t,x) = (iTThyf] [ [ [ ei(^(*'^''''^^-<2^'^>);iJ6j(t,x,?7,^) d9 dr] <p{y) dy + 0{W^) 
j=o -^^^ -^^^ 

(1.4) 

Here: 

k>C sup sup \d^d^n{x,p)\'^ (1.5) 

|q!|+|/3|>2 (x,p)elR2" 

for some C > 0. Moreover the following assertions hold: 

(1) S generates the graph At of the Hamiltonian flow 0^ : T*IR" T*W- \ft £ [0,T]: 

At := {{y,r];x,p)€ x T^M" \ {x,p) = 4>'^{y,r])} 

= {{y,7];x,p)eT*W'xn''\p = V^,S,y = Vr,S,0 = VeS} (1.6) 

(2) S G C°°([0,T] X M^" X R^;R) and has the expression: 

S = {x, r?) - - t{Lx, x) + {Q{t)e, 9) + {v{t, x, rj),e + f{t, x, 6)) + {iy{t, x, r/, 9), 6) 

+ g{t,x,v,9). (1.7) 

Here f{t, x, 9) : [0, T] x x R*^ ^ R'^, u{t, x, r], 9) : [0, T] x R" x R" x R'^ ^ R'^, g{t, x, rj, 9) : 
[0, T] X M" X M" X ^ R and Cafia{T) > are such that 

sup [\d^dU{t,x,v,9)\ + |5^9^9,^5(t,x, 77,^)1 + \d^d^^d^i^{t,x,v,e)\] < Co^^AT). 

[0,T]xR2n + fc 

The function (x,r?) v{t,x,rj) : R" x M" R is linear Vt G R, and t Q{t) : [0,r] ^ 
GL(A;) with Q(0) = 0. 

f^Sj The transport coefficients bj : j = 0, . . . are determined by the first order PDE: 
dtbo + ^V^5 V^bo + 2^ A^5 6o(i, a;, r/, 9) = Qn, 

j = (1.8) 

bo{0,x,v,9)=p{9). 

' dtbj + V^.6j + bj - -^A^6j_i = 0, 

2m 2m ^. > ^^_g^ 

^ 6,(O,x,r/,0) =0, p(0) G5(M^M+), = 1. 
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Here @n ^ C^(M + ;M) is arbitrary within the requirement: 

H^Gtv e C^{R^''+'';R), < a < N; 

(3)\/0<t<T,0<T< +00, the expansion ^1-4^ generates an parametrix of the 
propagator U{t) = e*-'^*/'^.- each term is a continuos FIO on 5(M") denoted Bj{t), j = 0,1, . . ., 
which admits a continuous extension to and: 

00 

^iHt/h ^ J2 Bj{t) + 0(^1°^). (1.10) 
j=o 

The notation 0{h°^) means: 

N 

\\RNmL^^L^<CN{T)n^^\\fN>0, VtG [0,r], ii^(t) :=U{t) -Y,Bjit). 

j=0 

Moreover, the expansion il.lO\) does not depend on p provided \\p\\l'^ = 1. Namely, if pi 7^ P2-' 

N N 

Y.B,[p,m-Y,B,[p2m = o{h''^')- 

j=0 j=0 

By applying the stationary phase theorem to the oscillatory integral ()1.4p . the integra- 
tion over the auxiliary parameters 9 can be eliminated and the WKB approximation to the 
evolution operator is recovered, necessarily multivalued on account of the caustics. 

Theorem 1.2. Let V{x) = ^Ix]"^ + Vo{x) with sup^g^n ||V^Vo(2;)|| < 1/ 'fh{il) be the h- 

Fourier transform of the initial datum ip. Then Vt € [0, T], t ^ (2r + 1)^, r G N, there exists 

Af 

a finite open partition M" x M" = Z)^ such that the solution of can be represented as: 

1=1 

ij{t, x)= [ Unit, X, ri) ^niv) drj, 0<t<T, t / (2r + 1)^ 
Un{t,x,rj) =Y^ e^^"(*'-'^) |detV^5(t, x, ry, 0* (t, x, 77))|-5eX-«6«,o(i, x, r/) + 0(n) (1.11) 



a=l 



Sa := S{t,x,r],ei^{t,x,7])), baft := bo{t,x,ri,9^{t,x,r])), aa := sgnVlSit,x,r],0^it,x,r])) 
where M is a t— dependent natural and: 

(i) On each Dn the equation = VeS{t, x, rj, 9) has £ smooth solutions 9^{t, x,r]), 1 < a < i. 
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(ii) Any function Sa{t,x,r]) solves locally the Hamilton- Jacobi equation: 

V x'-'a] 



2m 



-{t,x,r]) + V{x) + dtSa{t,x,r]) = 



(Hi) An explicit upper bound on the t-dependent natural M is computed in 12. 68\) . 

Example In the harmonic oscihator case V{x) = ^x"^ and the phase function is exactly 
quadratic S{t,x,r],9) = {x,r]) — ^{rj^ + x^) + {v(t,x,r]),9) + {Q{t)9,9) It admits a unique 
smooth global critical point 9*{t,x,r]) on {x,r]) £ M^" for t £ [0,T], t / (2t + l)f , r G N. 



Hence the series (jl.lip reduces to just one term conciding with the well known Mehler formula: 



ip{t,x) 



cos(t) 



Remarks 



1. The phase function is constructed (Section 2) through the Amann-Conley-Zehnder re- 
duction technique of the action functional ( |AZj . |CZj . [Car] ). Namely: 



S{t,x,r,,9) = {x,v)+ / hP{s)r{s)-H{j^s),r, + ^P{s))] ds 
Jo 



(1.12) 



where the curves T{t,x,9) = {'y^{t,x,9){s),^^{t,x,9){s)) are parametrized as follows: 



T{t,x,9) := < 



^-{t,x,9)is)=x- c^-it,x,9)iT)dT, cp- = 9''i-) + nt,x,9)i-), 



jP{t,x,9){s)= (tf{t,x,9){T) dr, cf = 9P{-) + fP{t,x,9){-) 
Jo 



(1.13) 

Here 9 G Pm-^^^([0, T]; M^") ~ M*-' (Fm is the finite dimensional Fourier orthogo- 
nal projector, k = 2n{2M -|- 1)) so that the parameters 9 can be identified with 
the finite Fourier components of the derivatives of the curves 7. ()1.12p represents a 
global generating function if k fulfills the lower bound (jl.Sp . In turn, the functions 
ir, F) ■■ [0, T] X X Pj\/L2 ^ Qa/L2 X Qm^^ are determined by a fixed point func- 
tional equation, essentially the Qm projection of the Hamilton equations (Section 2.3). 
The parametrization ()1.6p entails that S" is a smooth solution of the problem: 

( IV 

' ^-{t,x,ri,9) + V{x) + dtS{t,x,7],9) = 0, 

(1.14) 



2m 



S{0,x,r^,9) = {x,7]); Ve5(t, x, ry, 0) = 0. 
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2. Any function S{t,x,r],9) solving ()1.14p . i.e. the Hamilton-Jacobi equation under the 
stationarity constraint VgS = 0, is the central object to determine the so called geo- 
metrical solutions of the Hamilton-Jacobi equation (see for example the recent works 
|Carj . |B-Cj ). Global generating functions are clearly not unique and this is due to the 
presence of the 0- auxiliary parameters. Uniqueness holds instead for the geometry of 
set of critical points: 

S5 := {(x, r], e) e m2"+'= I VeS{t, x, t], 6) = 0} 

which does not depend on S because it is globally diffeomorphic to A^; a detailed 
study of Y^s is done in Section 2. We prove (Section 3) that symbols coinciding on S5 
generate semiclassical Fourier Integral Operators differing only by terms 0(/i°°). This 
will allow us to select symbols in such a way to make essentially trivial the proof of the 
continuity of the associated operator. 



3. The symbol 60 solving the geometrical version (jl.Sp of the transport equation is 

6o(t,x,r?,e) =exp|-^ [' A,S{T,-f^t,x,e){T),r,,e)dT] p{e) (1.15) 







If T2 > Ti, then k(T2) > k{Ti) so that r(ri,x,6') C T{T2,x,e). In the limit T ^ 00, 
6 ^ (j) £ L^(R+;M-^"') and we get the simplified functional (still well defined): 

bo{t,x,cP) = j\,V{x - jy^'iX) dX)clT^ p{cl)) (1.16) 

This corresponds to the zero-th order symbol of the Laptev-Sigal construction [LS]: 

Mt^y^v) =^^p|2^_^ ^xV{x''{y,r]))dT^ . 

Namely, the functional is the same, but is evaluated on the classical curves (with initial 
conditions x^{y,7]) = y, p^iy^r]) = if) instead of all the free curves used in ()1.16p . with 
regularity and boundary condition 7^(t,x,0)(t) = x. 

4. For potentials in the class ()1.2p and < t < T small enough no caustics develop, and 
there is a unique smooth solution 6*(t,x,r]) for (x,r/) G M?"^. The stationary phase 
theorem yields the 0-th order approximation to the integral (jl.4p : 

Uj:^\t, X, 7]) = ei^(*'^'^'^*)|detV^5(t, x, 77, r)|-5eX-6o(t, x, rj, 9^) + 0{h) (1.17) 

which coincide with the WKB semiclassical approximation. This fact suggests a rela- 
tionship, at any order in h, between the present construction and those of Chazarain 
|Chj and Fujiwara [Fu]. This is the contents of Theorem 14.21 



7 



5. The first three assertions of Theorem 1.2 represent the counterpart (in the rj variables) 
of a result of Fujiwara |Fuj . valid under the additional assumption that the number of 
classical curves connecting boundary data is finite. 

We thank Johannes Sjostrand for suggesting us the formulation of Theorem 1.2, and Kenji Yajima 
for a critical reading of a first draft of this paper. 



2 Generating functions for the graph of the Hamiltonian flow 

2.1 Lagrangian submanifolds and global generating functions 

Adopting standard notations and terminology (see e.g. [Wej ). we denote hy oj = dp A dx = 

n 

dpi A dx^ the 2~form on T*R" that defines its natural symplectic structure. As usual, 

i=l 

a diffeomorphism C : T*W^ — t- T*]R" is a canonical transformation if the pull back of the 
symplectic form is preserved, C*uj = lo. 

We say that L C T*W^ is a Lagrangian submanifold if uj\l = and dimL = n = ^dim T*M". 
In a natural way, a symplectic structure u) on r*M" x T*R" = T* (M" x M") is the twofold 
pull-back of the standard symplectic 2-form on T*M" defined as a) := prjo; — pr^uj = dp2 A 
dx2-dpiAdxi. Similarly, A C T*WxT*W is called a Lagrangian submanifold of r*R"xr*M" 
if (D|a = and dim(A) = 2n. 

A Hamiltonian is a C^-function H : T*W M and its fiow is the one-parameter group 
of canonical transformations 0^ : U C r*M"' — )• r*M" solving Hamilton's equations 7 = 
JVH{'^) {J the unit symplectic matrix) with initial conditions 7(0) = (xo,Po) ^ U . 
The Hamilton-Helmholtz functional: 

^[(7^7^)] := f[l''{sW{s)-n{^''{s),l^m ds (2.1) 
Jo 

is well defined and continuous on the path space i?^([0,t];r*M"). The action functional: 

^[7^:= f C{^'{s),r{s))ds (2.2) 
Jo 

is defined on H^([0,t];W^). In this paper we consider Ti = h V(x), so that the Legendre 

2m 

transform guarantees the corrispondence of the stationary curves of these two functionals. 

Definition 2.1. A global generating function for a Lagrangian submanifold L C T*M" is a 
function 5 : M" x M such that 
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o L = {{x,p) G r*R"| p = V^Six,e), = VeSix,e) }, 
o rank {Vl,S Vl,S) 



= max. 

L 



= max. 

A 



Similarly, a global generating function for a Lagrangian submanifold A C T*M" x T*M" is a 
C2 map S" : X E" X ^ M such that 

o A = {{x,p;y,r]) gT^'W xT*W\ p = V,S{x,r],e), y = V^S{x,r],e), = VeS}, 

o rank {vl^S V^.S Vl,s] 

It is important to remark that the foUowing set: 

Ss := {(x, r/, 0) G M" X M" X I = VeS{x, 9)] (2.3) 

is a submanifold of M^""*"'^' and it is diffeomorphic to A. 

We focus our attention on the graphs of a Hamiltonian flow 0^ : T^M" — )• T*M", which 
correspond to a family of Lagrangian submanifolds in r*IR" x T*M": 

At := { (y , 7?; X, p) G ™" X \ {x,p) = 4>\^{y,ri)] 

An important object in what follows is the family of global generating functions: 

kt = {{y,'n-x,p)Gn^xn^\p = V^S, y = V^S, O = Ve5(t,x,r/,0)} 

to be explicitly constructed for arbitrarily large times in the next Section. As is known, this 
technical tool has been developed in the framework of symplectic geometry and variational 
analysis (see |AZj . |CZ) . [Cha| . [LSik| . |Vit] . |Sik86| . |Sikj ) to sidestep the locality in time 
generated by the occurrence of caustics. 

2.2 Generating function with infinitely many parameters 

In the following we mainly review the construction of a generating function with infinitely 
many parameters described in [Car]. We begin by the following simple result (see |We| ) : 

Lemma 2.2. Let us consider the Hamilton- Helmholtz functional A[-] as in A curve 

7 g r(°) := {7 G H'^{[0,t];T*W^) \ 7^(0) = 0, 7^(t) = x} satisfies Hamilton's equations with 
boundary conditions: 

j = JVn{j), 7^(0) =0, 7"(t) = x 
if and only if the following stationarity condition of variational type holds: 

:^(7)[t;]=0 Vt;Grr(°) 
iJ7 



Proof. By computing the Gateaux derivative of the functional, we get: 
DA /"* 

— {j)[v] = / [7 - JVn{j)]{s)v{s) ds + Y{s)v''{s)\l yv G TrW. 
^7 Jo 

Now use the boundary condition 7^(0) = and recall that for TT^^^ it must be v^{t) = 0. 
The result is proved. □ 

The above Lemma has an important consequence: it allows us to introduce the notion of 
generating function with infinitely many parameters. 
First of all, it is easy to observe that the set of curves 

j{t,x,^){s) := (x- tr{r)dT, r <f{r)dT\ = ((/.^ 0^) (2.4) 



\ Js Jo / 

gives a parametrization of the path space T^^^ introduced in the previous lemma, namely: 

rW(t,x,</.) := {j{t,xA)i-) I <P e L2([0,r];M2")} 

Second, we define the functional with infinitely many parameters specified by E ^^([0, T]; M^") 
in the following way: 

Definition 2.3. 

S{t,x,rj,4>):={x,v)+ fb''{s)ris)-H{j-{s),rj + jP{s))]ds (2.5) 
Remark 

Introducing the traslated curves C = {C^, C^) := (7^) ^+7^); it is easy to see that the functional 
()2.5p admits the equivalent representation 



5(t,x,r?,</)) = {CiO),v) + feis)Cis) - H{C{s),e{s)) ds 







where now C^(*) = ^ C^(0) = V- 



Let us now make our assumptions on the Hamiltonian Ti more precise: 

Definition 2.4. ^ ^ 

nx,p) = ^ + V(.x) = ^ + {Lx, x) + Vo{x) (2.6) 

where Vq E C7°°(M"), L E GL{n), and 

\d^Voix)\ < Co, 

This allows us to look more closely at the structure of the generating function: 
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Lemma 2.5. The functional S admits the representation: 

S{t,x,ri,(l)) = {x,ri) - -^rf' -t{Lx,x) + {R{t)(j),(j)) + {v{t,x,r]),4)) +a{t,x,4)). 

(2.7) 

Here v(t, x, rf) has a linear dependence with respect to (x, rj) variables, and a{t, x, (f)) is bounded 
with respect to {x,(p). 

Proof. It is easy to see that: 

ft / rs 



S = {x,ri) + (^J^ r{r)dT^ r{s)ds 

- ^ + '^(^)^^) +l(^X- jy^T)dT^ ,(^X- (l)^{T)dT^ds 

- (p''{r)dT^ ds 

= {x,ri)-^r)^-t{Lx,x) + {R{t)(l),(l)) + {v{t,x,ri),(l))+a{t,x,^) (2.8) 

where 

ct>) := <f{r)dTr{s) - ^ 0''(r)dr) ' - L r{r)dr j'^ r{r)d^ ds 

{v{t,x,ri),(t>) := (^-^J^ ^{T)dT + 2Lx (t)''{T)dT^ds 

a{t,x,(f>) := - j\^[x- jy'^{T)dr^ 
The boundedness of a is immediate: 



ds (2.9) 



sup \a{t,x,(l))\ < t sup |Vb(2;)| 



Finally, we consider the orthonormal basis of L^, ea{s) = -^ct'"*, a eZ and the cor- 
responding Fourier expansion 0(s) = ^ (/)^°^eQ.(s). This entails the identification (j){s) = 



{</'^"^}oez £ ^ under the usual norm |^| = X^^I^^^^P generated by the scalar product 

(V','/') =EaV'("V(")• 
□ 

Proposition 2.6. The graph of the Hamiltonian flow 

At := {{y,v;x,p) e r*M" x T*R- \ {x,p) = cf>\^{y,r,)} 
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is generated by S: 

PS 

'd4> 



At = \ (y, v, x,p) G X ™" I p = V,S, y = V^5, 



DS 

Proof. The first component of the stationarity equation = — — reads: 

Dcp 

for all € L^. This is satisfied if and only if 

ris) = ^(^r^ + l^<fir)dT^ (2.10) 

that is 

nt, X, cl>){s) = - (7? + Y{t, X, cl>){s)) (2.11) 
m 

On the other hand, the second equation reads: 

= ;^('^)[^"] = (/' ^''i^)dr^ ^"(s) + Vy - </)^(r)dr) v%T)dT ds 
for all G L^. Integrating by parts, we get 

= f (t^{s)ds fv'-'{T)dT- /V(S) r v''{T)dT + VV (x - ^ (P"" {T)dA f v"" {T)dT ds 
Jo Jo Jo Jo \ Js J Js 

{T)dT ] / v^{T)dT ds. 



This entails 



/ (jfis) I v''{T)dT + VV fx - I 0^ 
Jo Js \ Js 

^{s) = -VV (^x - 0"^(T)dT^ (2.12) 

that is equivalent to 

Yit, X, 0)(s) = -V,y(7"(t, X, 0)(s)). (2.13) 
By a simple computation and ()2.10p we get: 

V„5 = x-t^- — r (l)PiT)dTds = x-t^- t ( (/)^Ys) - —7^ ds 
m m Jq Jq m Jo V ^ J 

= x-t^- ^''(s)ds + t^ = x- t ^^(s)ds = 7^(t, X, 0)(O) = y 
m Jo m Jq 

Finally, by (j2.12p . we can complete the verification: 

rt / rt 



= ^~ j^^^ (t>''{r)d'^ ds 



r]+ (j)P{s)ds = r] + 'yP{t,x,(l)){t) =p 







□ 



The following statement is a direct consequence of the above result: 
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Proposition 2.7. The Hamilton- J acobi equation is solved on the stationarity points 
more precisely S is a smooth solution of the problem: 



dtS{t,x,r],<j)) + 



2m 



it,x,ri,(j)) + V{x) = 0, (t,x)GM+x 



PS 

'd4>' 



(2.14) 



S(0,x,r],(j)) = (x,r]), ^^(t,x,v,(J)) = 0. 
For the proof we refer to [Car] (sections 3 and 4). 

Remark 2.8. As we have seen in Proposition \2.b\ fix t £ [0,T] and define the map 
Gt : M^" X L2([0,T];E2") ^ L2([0, T]; M^") 

Gt(x,7?,r,0^) :=f^ + - f ctf{T)dT,-Vv(x- fr{r)dT)). 

\m m Jq \ Js J J 

Then, the fixed point equation on L^([0, T]; R^"); 

(j) = Gt{x,r],(j)) 
is equivalent to the stationarity equation 

DS 

On the other hand, the solution of this equation determines the curves 

C{t, X, ^){s) := (x - r <A"(r)dr, r? + T <AP(r)dr) <A = (0^ c/f) 



(2.15) 
(2.16) 

(2.17) 



(2.18) 

\ Js JO y 

solving the Hamilton's equations ( = JVif(C) with boundary conditions C^(i) = X; C^(0) = rj- 
The following result deals with some topological properties for the set of the solutions. Let: 



2v Aq 



(2/3 + 1), a = l,2,...,n; /3 G N 



X{x,r]) := vT+IxP + IttP 



(2.19) 
(2.20) 



where A„ : a = 1, 2 



, ^, . . . , 



are the eigenvalues of L+L^. Remark that ta,i3 are just the resonant 



times of the hamiltonian flow generated by Hq := h {Lx,x). 

2m 

Proposition 2.9. There are D{T) < +oo, K2{T) < +oo, Ki{t) < +oo such that the 
solutions of equation { 2.17\ ) fulfill the estimates: 



L2 > K2{T)Xix,r]) VtG]0,r], \x\^ + \rj\^ > D{Tf; 
L2 < Ki{t)Xix,r]) Vt^t„,;3, V(x,r?)eM2". 



(2.21) 
(2.22) 
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Moreover, there is E{t) < +00 such that the difference of any two solutions (p, ijj of \2.11^ 
fulfills the estimate 

U-'^\\L2<E{t) yt^ta,^, V(x,r/)GM2". (2.23) 
Proof. We begin by remarking that the equation ()2.17p 

can be rewritten as 

- C{t, 4,) = ^o(i, X, ri) + X, (/)) (2.24) 

where 

£(i,</.) := (^^l^^ir)dT,{L + L^)l\-{T)dT^ 
^o{t,x,r]) := (^,-{L + L^)x) 

^i{t,x,ct)) := (^0,-VVo (^x - (l)''{T)dT^^ (2.25) 

To prove the inequahty ()2.2ip . remark that the non-degeneracy of L + entails the lower 
and upper bounds: 



Wo{T)X{x,ri) < \\^o{t,x,ri)\\L2=T-2 ( \\ + \{L + L^)x\^ ]' < Co{T)X{x,r,). (2.26) 



2 



Here Wo{T) := T2^j^^, Cq{T) := T'i;^^ and HM^l^m are the maximum and the minimum 
eigenvalue of the matrix: 

X = 



tI 



(L + Lt)2 



respectively. Moreover, 



\\^i{t,x,ct>)\\L2 = (^j^ \VVo(^x- j%%T)dT^ l^d^y <r^||Vyo||co =:C'i(r) (2.27) 
Now set M.{t, (p) := (p — C{t, (p). Hence the solutions of the equation ()2.24p fulfill the estimate: 

sup \\M{t,-)\\L2^L44>\\L^ > \\M{t,(j,)\\L2 = \\^o{t,X,r,)+^i{t,X,(l))\\L2 

te[o,T] 

> Wo{T)Xix,r])-CiiT) (2.28) 
For |xp + > D(T)2 we have Wo{T)X{x,r]) - Ci(r) > ^°^^\ {x,r]), and this implies 

L2>( sup \\M{t,-)\\L^^L^] ^^A(x,??) =: K2{T)X{x,r^) 
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Now consider equation ()2.24p in the particular case of Vq = (so that the Hamiltonian is 

p2 

T-Lq := h {Lx,x)). It becomes: 

2m 

(p- C{t,(j)) = ^o{t,x,r]) (x, r?) G X M'^ (2.29) 

The exphcit representation of the flow for the harmonic oscillator is 0^^, (Co ' Co) ~ ^''^(Co ' Co) 
where 



U 



I 



-(L + L-^) 

It is easy to prove that outside the resonant times ta^^ the flow can be globally inverted 
with respect to the boundary conditions {x,r]), namely: C^{s) = C^{t,x,r]){s) and C'^'i^) = 
C^{t, X, rj){s). By recalling Remark l2.8t this fact is equivalent to the existence of a unique global 
smooth solution (j)Q{t,x,r]) for equation (|2.29p . This argument works y{x,ri) G x M". In 
the particular case x = rj = ()2.29p reduces to: 

7W(t,(/)) :=0-£(t,0) =0 (2.30) 

The uniqueness of the solution implies that the linear operator Ai{t,-) : L^([0, t]; M'^") — )■ 
L^([0, t]; M^"), t 7^ ta,i3, is invertible. Now, we can come back to the general equation (j2.24p . 
written in the equivalent form: 

= M-\t, ^o{t, x, 7]) + ^i{t, X, </))) (2.31) 

for all t / ta,(5- By ()2.26p and (j2.27p . we get the inequality 

< \\M-\t,-)\\L2^M{T)X{x,ri) +Ci{T)) < Ki{t)X{x,r]) (2.32) 



where Ki{t) := \\Ai^^{t, ■)\\l^^l'^{Co{T) + Ci{T)). Finally, we have to prove the bound for 
the difference of any two solutions (j), if) for the equation (|2.3ip . In order to do this, we rewrite 
it under the form 

(t> - M~\t, ^i{t, X, 0)) = M-\t, ^o(t, X, 7j)) 

As a consequence, 

cj) - ^1; = M-\t,^i{t,x,^)) - M-Ht,^i{t,x,iP)). 
Recalling ()2.27p we have 

U-iph^ < \\M-\t,)h2^L^\\^i{t,x,(P)--f,{t,x,ij)\\L2 < \\M-\t,)h2^L^2Ci{T) =: E{t) 
and this concludes the proof. □ 
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2.3 Generating function with finitely many parameters 

In this section we describe how the global parametrization of the graph Aj of the Hamiltonian 
flow can be actually obtained through a generating function with finitely many parameters. To 
this end we use the reduction of the Hamilton-Helmholtz functional due to Amann, Conley 
and Zehnder (see jAZj . |CZj . |Carj . |B-Cj ). In this way we find, for the graph Aj of the 
Hamiltonian flow, a global parametrization of type: 

At = {{y,mx,p)en''xn^\p = V^S, y = V^S, = VeS{t,x,r],9)} 

The essence of the Amann-Conley-Zehnder reduction is the existence of an underlying finite 
dimensional structure for the equation investigated in the previous section: 

(0^ <t>P) = Gt{x, ri, ct>\cl>P) (0^ cjf) e L\[Q, T]; M2«). (2.33) 

We can indeed consider the two orthogonal projectors 

\r\<M \r\>M 

generated by any orthonormal basis of L^([0, T]; M^"); for instance er{s) := -^e~*^^ : r G Z. 
Then, let us introduce the decomposition: 

{no),ne)) = qMGt{x,7^,e- + r{e),ep + r{9)) (2.34) 
{9^en = FMGt{x,r,,e^ + n9),ep + ne)) (2.35) 

and prove the following 

Lemma 2.10. For M G N large enough the functional equation Iji2.34\ ) admits a unique 
solution f{e) : ¥mL^ ^ QmL^- The solutions of IQgl) can then be written in the form 

(0-,0p) = (r + r(0),e*' + r(0)) 

where 9 € Pm^^ — 0,1^^ finite dimensional parameters solving the fixed point equation 
^MW on IR^ k = 2n{2M + 1). 

Proof. Let us first verify that, if M G N is large enough, the equation (|2.34p realizes in fact a 
contraction on C'^(L^, L^). Hence it admits a unique solution f{t, x, 9) = {f^{t, x, 9), f^{t, x, 9)). 
By (|2.15p . (|2.16p the two equations read: 

(r (t, X, 9){s),fP{t, X, 9)is)) = Qm rmr) dr, -VV - j' 9^{t) + T (0)(r)dr 
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{9^{s),ep{s)) = p,^ ( ^ + 1 ep{T) + /p(e)(T) dT, -vv (^x - ^* r (r) + n9){T)dT 



m m 



It is proved in [Car] (Lemma 6) that the contraction property holds if: 

r2 sup \V^nix,p)\ < 1 (2.36) 

By Definition 2.4 it follows that sup(a; p)g7^*iRn \V'^'H{x,p)\ < +qo and consequently given 
< T < oo we get the contraction property for the first equation choosing M{T) large 
enough. In general, the second equation have many solutions depending on the values of 
{t,x). □ 

In this finite dimensional setting, we can consider the following set of curves, with t £ [0, T]: 

7"(t, X, e + f{9)){s) = x- 1' </<"(*, X, e){T) dr, ^"^{t, X, e) = e^ + r{t, x, e), 

< ' (2.37) 

7^(t, X, e + fie)){s) = r nt, X, e)iT) dr, cp^it, x, e) = ep + r{t, x, e) 



JO 

We note that this is a finite reduction of ()2.4p . but still contains all curves solving Hamilton's 
equations with boundary data 7^(t) = x and 7^(0) = because of (p are solving equation 
(1233]). Moreover, by a Sobolev's immersion theorem, T C i?H[0, T]; T*M") C C°([0, T]; T'^M") 
and this entails their continuity. 
We can now proceed to define the main object of this section: 

Definition 2.11. The finitely-many parameters generating function of At is defined as: 
S{t,x,rj,9) := {x , r]) + [\p (s)^ (s) - H {^^ (s) , v + Y (s)) ds 



■y{-)=-yit,x,e+f(t,x,e)){-) 

= Sit,x,r],e + fit,x,9)). (2.38) 

Here S is the infinite dimensional generating function of Definition 2.3. 

Remark The above generating function is fully parametrized hy 9 + f{t,x,9), 9 € M'^, and 
not by an arbitrary (p £ L^. This is the core of the finite reduction. 
Now we provide a more detailed study about the analytical properties of /. 

Lemma 2.12. Consider the pair of functions (f^ , f^) . Then: 
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(1) if^,f^) fulfill the following equations 



ris) = -Qm [ Qm [ {L + L^)f^r)drdT + <^>%t,x,e,f-)is) 
m Jo Jt 

Fis) = Qm [ iL + L^)f^iT)dT + <fP{t,x,e,f-){s) 

J s 

^^t,x,e,f^){s) := -^Qm j\MyVo{x- j\-{r) + f-{r)dr^dT. 
^P{t,x,e,f^){s) := -QmVFo (x-£r(r) + r(r)drj (2.39) 

(2) Under the condition 

d:=—\\QMf\\L + L^\\<l, (2.40) 
m 

they fulfill the estimates: 

\\f^{t,x,e){-)h2 < (i-d)-'^llQMfllvFollco 

\\fnt,x,e){-)h2 < T||L||||QM||||r(t,x,^)(-)||L^ + — IIQmIIIIVFoIIco (2.41) 

m 

(3) If in addition 

— IIQA/f sup sup|5^5^H(x,p)|2 < 1 (2.42) 

Nl + |j|>2 

then there exist Caa(T) > such that: 

\\d^dU{t,x,e){-)\y<CUT) (2.43) 

Proof. By direct computation, the first functional equation reads: 

r(s) = -Qm f fPmr)dT 
m 







\flM 1^ -QmVV (^x - 1^ 6^{r) + f^{r)dr^ dr 
\flM 1^ Qm{L + Lt) (^x - I' e%r) + f^{r)dr^ dr 
QAf VFo (x - f e%r) + r{r)dr\ dr 



m 

1 

m 







= -Qm Tqa/ f\L + L^)f^r)drdT + <!>'-{t,x,ej^){s) 
"^ Jo Jr 

wliere tlie last equality follows by ()2.39p . Analogous computation for f^i^s): 
fP{t, X, e){s) = Qm J\l + L^)f^{T)dT + cDf (t, x, 6, f-){s) 
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This proves Assertion (1). 

To see Assertion (2), remark that (|2.39p also entails: 



\^^{t,x,9,n{-)h2 < — IIQmIP llVyolbo (2.44) 

m 



whence we obtain: 



\\nt,x,e){.)h2 < — ||QMf||L + L^||||r(t,x,0)(-)k2 + ||cD(t,x,ry,0)(.)||i2 

m 

If we choose M large enough, then d := ^WQmW^WL + < 1 and hence we get 

\\nt,x,e)i-)h2 < {l-d)-'—\\QMf\\VVo\\co 

m 

In the same way we have the estimate: 

\\f{t,x,e){-)U2 < r^||L + Lt||||QM||||r(t,x,0)(-)k^ + ||<i>^(t,x,r?,e)(.)||i2 

< r5||L + Lt||||QM||||r(t,x,0)(-)k2+ri||QM||||vyolbo 

This proves Assertion (2). 

The equation for the first order partial derivatives reads: 

-j—{t,x,e) = ^ / Qm / {L + L^)^{t,x,e){r)drdT 
dxi m Jo dxi 



If M is large enough, then d' := gHQ^f ||i + ^^11 + SIIQa/ f llV^Folbo < 1 and we get: 

{t,x,d){-) 



ri 

L2 



dx. 

The equation for the second order partial derivatives reads: 
"J ^t,x,e) = ^ Qm / {L + L^)^^{t,x,e){r)drdr 



dxidxj ' ' ITT' Jo Jt dxidxj 



ni 



+ j\MFapk{t,x,e){T)(^5kj + ^£-{t,x,e){r)dr^ {^^'^ I ^-(*' ^) (0^^^^) 



where 



19 



As before, if we require d" := ^||Qm|P(||-^^ + -^^"^11 + l|V^Vb||co + \\F\\co) < 1 then 



,2 fx 



dxidxj 



<(l-d 



L2 



m 



F\\co {T-2 +T^Vn\i, +2T'\\Vrh2 



For the higher order derivatives in x and also for ^-partial derivatives we can proceed in the 
same way, with the general condition 

— IIQjuf sup sup\dldm{x,p)f <1 (2.45) 
in order to conclude the existence of Cacr{T) > such that 

\\d^d^fit,x,e){.)\\L2<CUT)- (2.46) 

This proves Assertion (3) and thus concludes the proof of the Lemma. □ 

Theorem 2.13. The generating function \2. 3^) admits the following representation: 



S 



{x, r?) - — r/2 _ t{Lx, x) + {Q{t)e, 9) + {v{t, x, r,),9 + f{t, x, 9)) + {u{t, x, 9), i 
2m 



+ g{t,x,9). (2.47) 
Here: 9 G R'^; t Q{t) G GL{n), Q(0) = 0; more overthere are Ca(3a{T) > such that 

The function {t, x, rj) i— )■ u(t, x, rf) is linear in x, rj, and finally: 



k > CT^ sup sup \d'^d^n{x,p)\'^. 

|a| + |/3|>2 x,p 

T nr 

Proof. By the explicit upper bound ||QAf|| < — y -77 of ()2.42p we have 

27r V M 

k>CT^ sup sup |a°5^^(x,p)|2. 

|a| + |/3|>2 x,p 

We recall the structure of the infinite dimensional generating function: 

5(t, X, ??, (/)) = (x, rj) - -^r?^ - t{Lx, x) + {R{t)(t), 4>) + {v{t, x, r]), 4)) + a{t, x, </>) 
Zm 

As a consequence: 

S{t,x,r],9) := S{t,x,ri,9 + f{t,x,9)) 

= {x, V) - i^V^ - t{Lx, x) + {R{t){9 + /(t, X, 0)), + /(t, X, 9)) 
+ (^;(^, X, r]),9 + /(t, x, e)) + a{t, x,e + f{t, x, 9)). 



(2.48) 
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We can thus make the identifications: 

{Q{t)9,e) := {R{t)e,9) 
{iy{t,x,e),9) := {2R{t)f{t,x,e),e) 

g{t,x,e) := a{t,x,e + f{t,x,e)) + {R{t)f{t,x,e),f{t,x,e)) (2.49) 

Now it is easy to see that 

\u{t,x,e)\<2\\R{t)\\\\f{t,x,e){-)h2<C{T) 

and this entails boundedness with respect to the 0-variables. Moreover the same property 
holds true for the other term. We have indeed: 

\9{t,x,e)\ < \\a{t,x,<P)i-)\\co + \\R{t)\\\\fit,x,e){-)\\l, < C'{T) 

By using the above results we then get the existence of C'{T) such that \g{t,x,6) < C'(T). 
The estimates for the partial derivatives follow in the same way. □ 

Theorem 2.14. The graph of the Hamiltonian flow 

At := {{y,v;x,p) G r*M" x \ {x,p) = </.^(y,r?)} 

admits a global generating function with finitely many parameters: 

At = {{y,7];x,p) en"" xn^" \p = V^,S, y = V^S, = VeS{t,x,7],e)} 

Proof. By Proposition (|2.6p we can write: 

f DS 

At = I p = V,5, y = V^S, = — 

where S is the infinite-dimensional generating function of Definition 2.3. Now we remark that 
the finite-dimensional stationarity condition: 

= v,5(t,x,7?,r) 

is equivalent to the variational equation expressing the stationarity: 

DS 

0=—{t,x,r],(l)') 

Indeed, by Lemma [2.10l and [Car] (see Lemma 7), there is a bijective correspondence between 
the solutions of the two equations, (p'* = 6* -\- f(t, x, 9*). Moreover it is easy to prove that 
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This is true because of the definition S{t,x,r],9) := S{t,x,r],0 + f{t,x,6)), and the compu- 
tation 

DS 

V^S{t,x,r],e) = xS{t,x,r],(t))\^=gj^fi^t,x,e) + -jj^it,x,r],(p)\^=g+f(^t,x,e)Nxf{t,x,9)]. 

Evaluating both sides on the solutions 6* we get the relation. The same argument applies to 
V^5, and this concludes the proof. □ 

Theorem 2.15. The Hamilton- Jacobi equation is solved by the smooth function S{t,x,rj,6) 
on the stationary points S5 = {{x,r],6) G Ji'^^+k | \/gS{t,x,r],9) = 0}. More precisely: 

'^"'^''-(i,x,r/,0) + F(x)=O, 

(2.50) 



dtSit,x,r],e) + 



2m 



^ S{0,x,ri,e) = {x,r]); (x, ??, 6*) G Ss. 

Proof. We recall that, by Proposition l2.71 the Hamilton-Jacobi equation is solved by S{t, x, r], 1 

DS 

on the infinite dimensional stationary points defined by -—--(t,x,r],(f)*) = 0. 

Dcp 



IV 5P 

dtS{t, X, r?, <j)) + ' ^ ' (t, x, rj, 0) + V{x) = 0, (t, x) G M x M" 



2m 



DS 



(2.51) 



5(0, X, ri, (j)) = (x.rj), -jr-{t, x, rj, 

D(p 



0. 



On the other hand, we have 

^ xS\{t,x,-q,e*) = xS\ it ^x,ri, <))*)■, dtS\{t^x,-q,e*) = dtS\{t^x,ri,(j)*) 



(2.52) 



Indeed the first equality is proved in the previous theorem; whereas for the second one we 
observe: 

DS 

dtS{t,x,r],e) = dtS{t,x,ri,(j))\^=Qj^fit-xfi) + jj^{t,x,r],(j))\^=e^f^t^^^Q)[dtf{t,x,e)]. 
DS 



Since -^{t, x, r], cj)*) = the second equality in (I2.52p is proved. ()2.5ip and (I2.52p then yield 
Dcp 

the assertion. □ 



Theorem 2.16. Let S and S5 be as in Theorem \2.15[ Then there exists Qn G C^([0,T] x 
]^2n+fc. 0jv|ss = 0, such that an equivalent generating function Sn is given by the 

solution of the problem 

'^""^"^''-(t, X, rj, 9) + V{x) + dtSNit, X, 7], 9) = On, 

(2.53) 



2m 

SN{0,x,r],9) = (x,r/). 
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n(ejv) := dive ( e^^^ ) (2.54) 



Moreover, defining: 

|V,5iv| 
Sn enjoys the property: 

n^'(e7v) G C5°°([o,r] X m2"+^;M) vi<i<iv, n = i,2,... (2.55) 



Proof. We remember that S5 C R '^^ is a submanifold of dimension 2n thanks to the 
nondegeneracy condition 

rk{VlgS, V^gS, VleS)\j:s = max = k = 2n{2N + 1) 

for some TV > 1. We define z := (x, 7/, 6) £ M2"+'= and for any point z E S5. Define furthermore 
S (not necessarily uniquely) through the conditions: 

dtS{t,z) = dtS{t,z) + L{t,z), (2.56) 
Vj{t,z) = V,S(t,z) + F{t,z), (2.57) 
Aj{t,z) = A^Sit,z) + G{t,z), (2.58) 

where L = (L^,L^,L'') e C^{[0,T] x E2"+^';M) and L{t,z) = 0, the perturbation of the 
gradient in (f2371> is F = (F^,F'',F^) G Cj°°([0,T] x M^n+fc. i^2n+fc) ^-^j^ F(t,z) = while 
in ()2.63p we require G G C^{[0,T] X M2n+fc.]^) g^^^ Q^^^^^ ^ Q addition, we require that 
F^(t, z) ^ for z ^ S5. Hence the new stationarity equation: 

VeSit,z) = F''{t,z) (2.59) 

implies = Sg. In order to verify (j2.54p we require a suitable asymptotic behaviour of 
L,F,G around S5. Indeed, 

dtS{t,z) = dtS{t,z) + L{t,z), VJ{t,z) = V^S{t,z) + F%t,z). 

So, by easy computations and by ()2.50p . we have 

e(t,z) = \^^{t,z) + V{x) + dtS{t,z) 

= 7^1 Vx5(t, z) + F"(t, z)|2 + y(x) + dtS{t, z) + L(t, z) 
2m 

= 7^1 V.5(t, z)|2 + V{x) + dtS{t, z) + -V^S{t, z)F^{t, z) + -^|F"(t, z)|2 + L(t, z) 
2m m 2m 

= -V,.S(t,z)F-(t,z) + -^|F-(t,z)|2 + L(t,z) (2.60) 
m 2m 
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Now we can always require that the vanishing asymptotic behaviour of F''^,F^,L around 
are such that it holds: 

n(e) := div (^0^^) e C^{[0,T] X r2-+^m) 

By the same arguments as above, we can look for Sn such that T,Sj^ = S5 and 

dtSN{t,z) = dtS{t,z)+LN{t,z), (2.61) 
V,SNit,z) = V,S{t,z) + FN{t,z), (2.62) 
A,SNit,z) = A,Sit,z) + GNit,z), (2.63) 

where F^,F^ and Lj\[ are chosen in such a way that: 

n^'(e7v) G c^{[o,T] X M^n+fc.^) yi<j<N 

□ 

Let us examine the topology of the finite-dimensional critical points set. 

Theorem 2.17. Let S be as in Definition \2.11[ and X{x,r]), ta^p as in {2.19]) and h2.2U\) , 

respectively. Consider {t,x,r]) g]0,T] x M" x W\ Then: 

(1) All solutions G Pm-^^^([0, T]; M^") ~ M'' of the stationarity equation: 

= VeS{t,x,rj,e) 

fulfill the estimates 

\\9\\l2 > K2iT)X{x,r]) VtG]0,T], \x\'^ + \r]\'^ > D{Tf- (2.64) 
||^||l2 < Kiit)X{x,T]) yt^t^^fs, V(x,r?)GM2". (2.65) 

where D{T) < +00, K2{T) < +00 while Ki{t) < +00 is a constant defined for t 7^ ta,i3 ■ 

(2) The difference of any two solutions 0,uj fulfills the inequality 

\\e - uj\\l2 < E{t) Vt/t„,/3, V(x,r/)GM2n_ ^2.66) 

Proof. By Proposition 12.91 and Lemma 12.101 all solutions (p G L^oi the variational equation 

DS 

= —{t,x,r],(f)) 

are such that 

ci) = d + f{t,x,e) 
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and fulfill the inequalities ([MI]), ([2:22]) and ([2:23]) for some constants Ki{t),K2{T) and E{t). 
Using the uniform bound proved in Lemma 12.121 

\\fit,x,e){-)h2<Coo{T) 

we easily establish the existence of the new constants Ki{t), K2(T) and E{t). □ 

Theorem 2.18. Let us suppose V{x) = + Vq{x) with sup ||V^Vo(a;)|| < 1, {t,x,r]) G 
]0,r] X M^" and t / (2t + r G N. T/ien i/ie following quadratic form: 

{VlS{t,x,r],9)u,u), uGR'', 

is non degenerate on all points solving 'VgS{t,x,r],6) = 0. 

Proof. The quadratic form is non degenerate if and only if the solutions 6 of the equation 

VeS{t,x,rj,e) = (2.67) 

are isolated points. This property can be translated in the infinite dimensional setting of the 
equation 





thanks to the equivalence (j) = 9 + f{t, x, 6) shown in Lemma 12.101 We recall that 



JO 



5 = (x,r?)+ / (/ 4>^{T)dT,r{s)) 



^{^ + <^P(r)dr^ -v(^x- ^'={T)dT^ ds 



Now we perform the partial reduction of the infinite dimesional parameters, by means of the 

DS 

first stationarity equation -——(t,x,r),(f>) = corresponding to mcjf(s) = rj + (IfMdT 

Dct)'' Jo 
(essentially, the Legendre transform). Therefore we get the new functional: 



S{t,x,r^,cj)n = {x- / </,-(r)dr,r/) + 



m 



'■^(s)|2 - y ( X - / ^''{r)dr 



ds 



Setting 7^'(s) = x — (j)^{T)dT, we can consider the equivalent form 



^[7l = (7^(0),^)+ / - \r{s)\'-Vij^is)) 



ds 



with the boundary conditions: 7^(t) = x and m7^(0) = rj. The second variation is: 



D'^A 1 



m [|(57(s)P - VV(7^(s))<57(s)<57(s)] ds 



Writing down the integrand under the form 
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-W(7^(s)) 







ml 



5j{s) 



we realize that requiring V'^V{x) non-degenerate Vx G M", then the second variation is a 
bihnear non degenerate functional. This implies that all the stationary curves of the action 
functional, namely the curves solving 

DA. 



are isolated points belonging to H^{[0, t]; M"). We conclude that the same property must hold 
for the points S R*^ solving equation (|2.67p . □ 

Next, we investigate the number of solutions of the stationarity equation. 



Theorem 2.19. Let us suppose V{x) 



1 Ul2 



+ Vo{x) with sup ||V'^Vb(x)|| < 1, {t,x,r]) e 



]0,r] X M^n i ^ (2r + l)f , r G N. Then the stationarity equation 

VeS{t,x,v,0)=O 

has a finite number of solutions 9^{t,x,rj), 1 < a < J\f{t). The upper bound has the expression: 

{2E{t)Y 



M{t) < 



e{TY 



(2.68) 



Here E{t) as in Theorem 2.11 whereas 

1 



e{T) :- 



k 



inf (t,x,77,e) supij 




{t,x,r],6) 




^^P{t,x,r],e) SUPj j„j 


83 S 


{t,x,r],e) + 1 





(2.69) 



Proof. By Theorem ()2.17p all critical parameters must be contained in the compact set Br C 
R*' with r := 2E{t). As a consequence, there exists a subsequence {Oa[j)}j<m converging to 
some point 6 in i?j,(0). However the function V QS{t,x,rj,-) is continuous on M.^ . Hence the 
limit is also a critical point, namely = 'VgS{t,x,r],9). By the previous theorem all the 
critical points of S are isolated. This is a contradiction, so their number must be finite. In 
order to obtain an upper bound for this number, we first observe that 

vls{t,x,v,e) = vls{t,x,ri,9*) + ^vls{t,x,v,e* + x{9-e*)) dX 

= VjSit, X, 7], 9*) + r DeVlS{t, x, t], 9* + X{9 - r)) dX {9 - r) 
Jo 



26 



We know that, thanks to Theorem 12.181 the first matrix on the righthand side is non dege- 
nerate. In order to verify that the addition of the second one does not change this property, 
we estabhsh the matrix norm inequahty: 



/ DeVls{t,x,r],e* + \{e -e*)) dx {9-e* 

Jo 



< \\VlS{t,x,r],e* 



(2.70) 



Here || • ||2 is the usual norm for the matrix viewed as an operator. Now denote e := \\9 — ^*||. 
The above inequahty is a fortiori verified if: 



e Vk sup sup 

{t,x,ri,8) i,j,m 



{t, X, r], 9) < —= inf sup 

yjh {t,x,ri,e) ij 



dOidOj 



{t,x,rj,e) 



(2.71) 



because the l.h.s is an upper bound for the l.h.s. of (j2.7Up and the r.h.s a lower bound for the 
r.h.s. of (|2.70p . ()2.7ip is in turn a fortiori verified if: 



e Vk sup sup 

\{t,x,r},9) i,j,m 



(t, X, 77, 6*) + 1 I < —= inf sup 



dOidej 



it,x,7],e) (2.72) 



and this yields (j2.69p . In this way, we have found the radius e(T) of the balls in M'^, where 
each 6* is a unique local critical point. This local confinement of critical points together with 
the global one proved in Th 12.171 allows us to get an estimate of their total number A^. We 
simply compute the ratio between the volume of the ball Bj. containing all the points and 
the volume of the small isolating balls. 

vol{Br) {2E{t)f 



M{t) 



V0l{Be 



□ 



We use Theorem 12.191 in order to study the global behaviour of the stationarity equation. 

Theorem 2.20. Let us suppose V{x) = \\x\'^ + Vq{x) with sup ||V^Vo(x)|| < 1, {t,x,r]) G 
]0,T] X ]R2" and t / (2r + l)f , r G N. Let the number J\f{t) he given by ^MW- Then there 



N{t) 



exists a finite open partition M^" ~ U such that the equat 



ion 



= VeS{t,x,ri,e) 
admits on each Di exactly I smooth solutions 0^{t,x,ri), 1 < a < £ 



(2.73) 



Proof. We recah that T.s := {{x,r],9) G M^"+^| = V0S{t,x,rj,6)} is a 277-dimensional 
submanifold of K^n+fc diffeomorphic to A^. Moreover, by the nondegeneracy hypothesis on 
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V'^V we have the transversal behaviour of (with respect to (x, r/) G R^") almost every- 
where; namely the rank of VgS can differ from its maximum value (k) only on subsets whose 
projection on (x, rj) £ M^" is of zero measure. The condition of transversality is fulfilled on 
components (locally diffeomorphic to open sets of M^*^) where the local smooth inversion 
of equation ()2.73p is possible, yielding i functions 9^{t,x,ri). This argument works up to the 
finite maximum value A/'(t). □ 



2.4 Transport equations 

We conclude this section by introducing transport equations in a global geometrical setting. 

Theorem 2.21. Let us consider p G 5(M*'';R) with \\p\\l'^ = 1. The transport equation written 
on the stationary points T^s of the generating function S, 



dtbo + -V^S V^bo + ^A^S bo{t, X, r,, 0) = 0, 
m Zm 



(2.74) 



, bo(,0,x,ri,e) = p{e), (x,7?,e) 
admits the following solution: 

bolt, X, 7], 9) = exp j'^ A,S(t, 7"(t, x, e)(T), r?, 0)dr| p{e) (2.75) 

where 7^ is the family of curves defined in \2.3T^ . 
Proof. The inital condition is immediately verified: 

60(0, X, r?, 0) = /5(e) 



Recalling the results of Proposition 12.61 and Theorem 12.141 we compute the expression of 

the differential operator dtb^ -\ V^jS 'Vxbo(t,x,rj,6) when evaluated on the submanifold 

m 

S5 := {{x,r],e) £ M^n+fc | q = V9S{t,x,rj,9)}. Namely, 



dtbo + —VxS V^6o(i, x, rj, 6) 
m 



1 



dtbo + - (r/ + ^""{t, X, e){t)) V^6o(t, x, r], 6) 



m 



dtbo + Tit, x,9){t)Vxboit,x,7j,0) 

d^bo{p., X, r/, 6) + 7^(t, x, 9){p)V ^boifi, x, r?, 6) 
d 



dp 



bo{p,j^{t,x,e){p),r],e) 



fl=t 



fl=t 

2.76) 
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where the expression of AxS{t,x,ri,9) 6o(t,x,r/, 0) is: 

2m 

^A,S{t,x,7],9) bo{t,x,r],e) = ^AxS{t,j''{t,x,9){fi),rj,e) bQ{^i,J''it,x,e){pi),r],e) 
Now, we write down the equation in the new variable /i and for all {x,ri,9) G M?"'~^^: 



|^6o(/i, 7"(t, X, e){fi),r^, 6) + ^A^Sifi, j^t, x, 0){fi),T], 6) bo{fi, ^^t, x, ??, 9) = 

If we define a(^) := bo{fi,'y^{t,x,9){fi),7],0) we can rewrite the previous equation as 

d 1 

_a(/i) = -— A^5(/i,7^(i,a;,0)(/i),r/,^) a(/i). 
a/x 2m 

where the variables {t, x, rj, 9) have to be considered as fixed. This yields: 

6o(//, 7"(t, ^, 0){^,),r,, 9) = exp 1--^ T A,S(r, 7^(t, rr, 0)(r), r?, 0)^x1 pi9) 







2m 

Finally, we make fj, = t and so we obtain the solution of the original problem (|2.74|) : 



1 



6o(t,x,7?,e)=exp<!-— AxS{T,j^t,x,9){T),r,,9)dT} p{9) 



□ 



Theorem 2.22. Let bo be defined as in [K7^ ) with p{9) := e~\'^\'' i{9) and i € C^°°(M^M+). 
Then boit,x,r],9) G C°°{[0,T] x M2n+fc.]^+-) bo{t,x,r],-) G 5(M'=;M+) for every {t,x,'n) 
fixed. Moreover, there exists a constant C^(T) > such that 

\d^boit,x,r],9)\ < C7+(r)e'^"™^''')e-l^l' V(x, r/, G R^n+fc ^2.77) 
Proof. Let us first obtain a more explicit expression for {Ax)S{-): 

A.xS{t,x,7],9) = 2tr{L)t +{Axu{t,x,9),9) + {v{t,x,r]),Axf{t,x,9)) 

+ 2{VMt, X, v),^xf{t, X, 9)) + A^git, x, 9) 

= 2tr{L)t + {2R{t)A,f{t, x, 9), 9) + {vH, x, r/), A,/(t, x, e)) 

+ 2{VMt, X, v), Vxf{t, X, 9)) + A^git, x, 9) 

= 2tr{L)t + {2R{t)A,fit, x, 9), 9) + {vH, x, r?), A,f{t, x, e)) 

+ 2{VMt,x,ri),Vxf{t,x,9)) + Axg{t,x,9) (2.78) 



Now, we recall that 



ft 

7^(t,x,0)(T) =x - / rXr) +r(t,x,r)(r) dr 
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where / and all its derivatives are uniformly bounded, as proved in Lemma 12.121 whereas 
V is linear in {x,r]) and g is L°° bounded. Now we observe that by setting p{9) := e~^^^^^{9) 
with a bounded ^ G C°°(M'^;M+) then bo{t,x,r], •) is a Schwartz function on M'^. Indeed, 

|&o(t,x,r?,0)| <exp|-L r|Ax-5(T,7"(t,x,e)(T),r/,0)|dr|e-l^l'eW 
But by the above detailed computation we see that 

\A,s{T,^^t,x,e),T],e)\ 

< |2tr(L)t| + 2||i?(t)|| ||A,/(^,7^e)||^2||0|| + ||7;(^,7^r?)||L2||A./(^,7^^)||L2 
+ 2||V,t;(t,7",r/)||i2||V,/(t,7",^)||L2 + \\A,g{t,j^ ,9)\\l^ 

\\AJ{t,^-{t,x,9),e)\\L2\\9\\ is linear in 9 and \\v{t,^-it,x,9),r])\\L2\\AJ{t,j-{t,x,9),9)\\L2 
has a linear uniform growth on (x, 0). To see this, remark that 

\\v{t,j^t,x,9),r,)U2 < ||7;(t,x,r/)||i2 + ||^;(t,£r(r)dr,r?)||i2 + ||z;(t,^V"(t,x,e")(r)dr,,7)|U^ 

The first and third term on the right hand side generate a linear growth on (x, ?]), the second 
term has a linear dependence on 9. The other terms above are bounded with respect to all 
variables. We conclude that bo has a uniform exponential behaviour on all its variables, and 
that the function p{9) := e"'^'^ ||^||(70 makes the effective dependence on 9 of Schwartz type. 
□ 

Theorem 2.23. Let S and Ss be as in Theorem\KT^ Then there exists Qn G C^{[0,T] x 
M^"-+'^;R) with @n\es ~ ^' ■^^^c/i that the solution Sn of the problem 



dtbo.N + —^xSn ^xbo,N + -^A^Sn bo,N{t, x, r], 9) = Qn, 
m 2m 



(2.79) 



bo,NiO,X,7],9)= pi9). 
enjoys the property: 

n^'(e^) G Cfe°°([o,r] X m2"+'^;M) vi<i<iv, iv = i,2,... (2.80) 

where, as in Theorem \2.16i 



n(Gjv) := dive Gjv7„ , 



Proof. Let us define 

5o,^(t, X, r?, 9) := exp A,Sn{t, ^^t, x, 9){t),v, 9)dT^ p{9) (2.81) 
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and prove that is solves the above problem. Indeed, we can write down the expantions for 
z := (x, ry, Q) around z G = S5 

6o,7v(t,2) = 6o(t,^) + /jv(t,^) (2.82) 
V^6o,7v(i,^) = V^6o(t,^) + 5Jv(t,2) (2.83) 
dtbo,N{t,z) = dMt,z) + hN{t,z) (2.84) 

all these terms are related to the choice of Gjv in Theorem 12.161 and their rate of convergence 
to zero near z are related as well. By unperturbed equation ()2.74p . we compute: 

0Ar(t, z) = dtbo,N + —VxSn ^xbo,N + -^A^Sn bo,N{t, z) 

m 2m 

= dMt, z) + hN{t, z) + — (V^5(t, z) + z)) {V^bo{t, z) + g^{t, z)) 

m 

+ -^-(A^5(t, z) + Gjv(t, z)){bo{t, z) + /^(t, z)) 
Zm 

= hN{t, z) + —VxS{t, z)gN{t, z) + F^{t, z) (V^6o(t, z) + qnH, z)) 
m 

+ l-A,S{t,z)fNit,z) + ^GNit,z){boit,z) + Mt,z)) (2.85) 

Moreover, by Theorem 12.161 we remember that around z G '^Sn ~ '^s the new stationarity 
equation is 

VeSN{t,z) = F^^{t,z) 

Now we can state that a suitable choice of F^,F^ and Gn leads to the following property: 
W{eN) G C^{[0,T] X M^^+'^jM) VI < i < iV 

□ 



3 A class of global FIO 

In this section we follow the general setting of Hormander |Hoj , and in particular we study a 
class of FIO related to the Hamiltonian flow 0^ , by using the generating functions constructed 
in the previous section. The study of the topology of their critical points will be useful to 
determine important analitical properties of the FIO such as asymptotic behaviour of the 
kernel and L^-continuity. 

3.1 Basic definition and main properties 

First, we introduce the set of phase functions: 
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Definition 3.1. The set of phase functions S{t,x,r],e) : [0, T] x M" x M" x M'' ^ M is the 
set of smooth global generating functions of the graphs At C T^M" x T*W^ of the canonical 
maps 0^ : r*M" — )• r*M", with the inititial condition S{0,x,r],6) = {x,r]). Each At admits 
the parametrization: 

At := { (y, r?; X, p) G T^R" x T*R^ | (x, p) = 0*^ (y , r?) } 

= {{y,r];x,p) Gn^" xTV \ p = V,S, y = VrjS, = VeS{t,x,r],e)} 

Before going further, we recall that by Theorem 12.171 the generating function S enjoys an 
important property. Namely, consider the set of critical points 

Ss := {{x, v, e) G m2"+^ I = VeSit, X, ??, 9)]. (3.1) 

Then is a manifold globally diffeomorphic to At; moreover for all t > the following set 

Ts:={(x,r?,0) GM2"+'= | \x\'' + \r,\^ > D{T)\ |0| < i?2(r)A(x, ry)} (3.2) 

is free from critical points, i.e.: 

Second, we introduce the relevant class of symbols associated to S: 

Definition 3.2. The set of symbols consists of all b G C°°([0,r] x M^" x R'';R) such that 
(^) 

6(0, X, 77, 6) = p{9), p{-) G 5(R^ R+), / p{e)de = 1. 



(a) For all multi-indices a,f3,a and t g]0, T] the inequalities 
\b{t,x,rj,9)\ < < 



^ C+(T) e^(^''')e-l^l', {x,r],0) ^ T5 



(3.3) 



C-(r)A-"(x,r/)e-l^l , (x,77,0)gTs 



hold for some constants g ^{T) > 0. 



'■,13,' 

Remark 3.3. The exponential upper bound outside Tg is verified by the symbol 60 (see 
Th. I2.22P and also, as we will see, by any other symbol bj,j = 1, . . . entering in Theorem 
1.1. Moreover, on domain there are no critical points for the function S and this leads 
to require asymptotic vanishing behaviour of type A~"(x, r/)e~l^l^ in this region for bo; as a 
conseguence the same asymptotic property is fulfilled by all bj. This setting is motivated by 
the fact that the contribution of this region to the FIO can be of order 0(/i°°) as we see in 
Corollary 13.71 In this framework, we provide a very simple proof of global continuity. 
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Finally, we introduce the class of global FIO associated to the Hamiltonian flow (pi^: 

Definition 3.4. Fix a phase function S as in Definition \3.1\ and a symbol b as in Definition 
Then the global h-Fourier Integral Operator on 5(M") is defined as: 

B{t)if{x) = {2Trhy f f f ei^'^^''^''''^'>-^y'i^'> b{t,x,rj,e) d9 dr] ^{y) dy (3.4) 
Jr" Jr" Jr'' 

In equivalent way, it can be rewritten in the form: 

B{tMx) = (27rn)-" [ [ efMt^''^y^^^ b{x,u) du ip{y) dy (3.5) 
Jr" Jr*" 

where u := {rj,6), S{t, x, y, u) := 5(t, x^ rj, 6) — {y, rj) and b{t, x, u) := b{t, x, r], 9). Indeed, if S 
generates the Lagrangian submanifold A, then S does the same in new variables: 

A = {{x,p;y,r]) £T*W xT*^' \ p = V^S, 7] = -VyS, = V„^} 

Theorem 3.5. Let us consider the FIO as in Definition \3.4\ 

B{t)ip{x) = {27Thr'' f f [ ei'^'^'^''''''''^^-^y'''^^ b{t,x,r,,d) de di] ^{y) dy 
Jr^ Jr" Jri^ 

Then B{t) : 5(M"') — )■ 5(M") is continuous and admits a continuous extension as an operator 
in L'^{W). 

Proof. We begin by rewriting the FIO under the form of an integral operator acting on the 
/i-Fourier transform of the initial datum: 

B{t)ip{x) = {Inhy [ [ eTi'^^^''=''^''^H{t,x,r],e) deipni7j)d7j 

Jr" Jr'' 

= (2^n)-" [ an{t,x,rj)(prMdv- (3-6) 

Mt,^,v) ■■= [ e^nSit.-.vfi) i^t^^^^^e) de 

JR* 

This is because of the integral in the ^-variables is absolutely convergent since b{t, x, rj, •) G 
5(M'^), and ip{y) is also a Schwartz function and therefore admits a /i-Fourier transform in 
5(M"'). The absolute convergence of the integral, as well as the L^-continuity, is the conse- 
quence of the following computations. 

Mt,x,7^) = I ei^^*'-'"'^) 6(t,x,r?,^) d^ = ei<-'^)[a+(t,x,7?) + a,7(t,x,,?)] (3.7) 
d^{t,x,ri) = [ e*(^(*'^'^'^)-<^'''>U(t,x,77,e) (3.8) 

a+{t,x,r]) = [ eE(^(*'^'^'^)-<^'^» 6(i,x,r/,0) (3.9) 
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with 6 := K2{T)X{x,r]). For t = we have a'^{0,x,rj) +?^(0,x,r/) = 1, B{Q)ip = 99, and the 
continuity is obvious. For t > we can apply the estimates of Property (ii) of Definition 13.21 
In the region containing the critical points we have: 

\dt{t,x,r])\ < [ \b{t,x,r],9)\de< [ C+{T)e^'^'''''^e-^^^'d9 

= Co+(r)e^(^''') / e-^^^'de (3.10) 



By writing down the integral in spherical coordinates, we have the following simple estimates 

f POO fOO 

/ e-\'\"de = cu e-P^p'-'dpiCcML) e'P'^dp = cML)e-''' 

Juf^XBsiO) Js Js 

for all L > and dk{L) := supp>Q e~ '^^ p^~^ e^^ . In particular we choose L := 1 + K2^{T), so 
that it follows 

\at{t,x,^)\ <C+{T)e''^''^^^cudu{L)e-^^^'^^^^^^^^ (3.11) 
Whereas in the other region we can write: 

|a,-(t,x,r?)| < / \h{t,x,r,,e)\de < I C^{T)\-'\x,r^)e-\'\"de 

Co"(r)A-"(x,r?)e-l^l'd0 = Cq {T)Tr-1 X-'^ix,!]) (3.12) 



< 



We can now apply the Schur Lemma to both integral operators and this yields the L^- 
boundedness. □ 

Now we prove a result, based on an argument of Duistermaat (see [Dui| . Prop 2.1.1). 

Lemma 3.6. Let us consider a FIO of type ^cl.4\) with phase function S and symbol g leading 



to a convergent integral. Suppose that 

Uig) := dWg (5^^) e C-([0,T] x M^n ^ ^k.^-^ ^g^^g) 

and that the FIO with symbol Il{g) is convergent. Then, the following equivalence holds: 
^lS{t,x,rj,9)g(^f^ X, r], 9) de^{r])dr] = -ih [ [ e:s^(*'^''?'^)ng(t, x, rj, 6) de0{r])dr] 

Jr" Jri^ 

(3.14) 

P™/. Th. diffe.ent.al ope.ato. :^ venfie. the .elat.on 

-ihhei^ = eTi^ 
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Now, by using integration by parts and the definition of the operator L, we get: 
I eTi^gcie = -ih I L (e^'^) gd9 = -ih [ e^^n(c/)d(9 



□ 



Corollary 3.7. Let g G C°°{[0,T] x R2n+fc.]^) ^y^/j f-f^^t U^(g) g C~([0,T] x M2"+'=;R) 
and that the corresponding FIO is convergent for all < j < N . Then: 

|5(i,x,„,e)~^^^ X, 7], 9) deip{-n)dr, = i-ih)^ [ [ ei^(*'^'''''^)n^5(t, X, r), 9) d90{r])dri 
Proof. The iterated apphcation of the previous Lemma gives the result. □ 

Remark 3.8. If we take two symbols 51, 52 coinciding on and moreover such that 51 — (72 = 
g, g as in the above Corollary, then the related FIO coincide up to order 0{h^). 

4 Global parametrices of the evolution operator 

Here we prove the main result of this paper. Consider the initial- value problem for Schrodinger 
equation 

ihdttpit, x) = -—A7p{t, x) + V{x)ip{t, x), 

2m (4;^) 

^^(0,3;) = (p{x) £ S{R''). 
with a potential V quadratic at infinity, of the type ()2.6p . 

We proceed to apply the results of the previous two sections in order to prove Theorem 1.1; 
namely, to construct a parametrix for the evolution operator under the form of series of a 
global FIO such that the solution of the Schrodinger equation (j4.ip admits the following 
representation: 

00 „ 

i;{t,x) = y{2TThy / e^(^(*'^'''''')-<J''''>)fi^6j(t,x,r/,^) d9 dr, <f{y) dy + 0{h^) 

within the time interval t £ [0, T] with T arbitrary large. 
Proof of Theorem 1.1 

Denoting ■= Ax + V(x) the action of the Schrodinger operator we look for a family of 

2m 

global FIO {Bj{t)}jfzfq with symbol bj{t,x,r],9) enjoying Properties (i) and (ii) of Definition 
31 such that 

00 „ 

= {H^- ihdt) Y / eE(^(*'^'''''')-<2^'''»;iJ6j(t,x,r?,^) d9 dr] ip{y) dy + Oih"^) 
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First of all, the approximation of order zero is the operator Bo{t) defined as: 

Bo{t)^ := {2TThy [ eE^^^^'^'^'^^-^J'-^^^boCt, x, r,, 6) dO di] ip{y) dy (4.2) 

JK2n + fc 

It has to reduce to the identity for t = and to represent the semiclassical approximation 
of the propagator. To this end, the related phase function solves the H-J problem ()2.53p and 
moreover the symbol bo solves the regularized geometric version of the transport equation as 
in Theorem 12.231 As we observed in Remark 13.31 ^6 require that in the region free from 
critical points of 5, the symbol 6o behaves as A~"'(x, ryje^'^l^ . 
Now, we easily see that 



{H^ - ihdt) ei^ho 



1 Q 



2m 



+ v{x) + dtS ] bo - in" ( dtho + v^bo + 



2m 



7;—^xbo 
zm 



The first two symbols of this sum vanish on the critical points set in such a way we can 
apply Corollarv l3.7l and so they realize bounded operators of order 0{h°°). As a consequence, 

e'i"^ bod&ipn{'n)d'n = — 

In + k 

The operator Bi{t) and the related symbol bi fulfill the analogous relation: 



{H^ - ihdt) [ ei'^bodB^rMdv = f ei^ A^bodeipn{v)dv + 0{ff^). (4.3) 



{H^ - ihdt) eT^^hbi 



2m 



+ v{x) + dtS bi - ih^ dtbi + v^s v^6i + 



2m 



2m 



Accbi 



The transport equation we now require for this symbol is the following: 

dtbi + Vxhi + 7^A^5 bi = ^A^bo, 
2m 2m 



(4.4) 



bi{O,x,r],0) =0. 
As a consequence, 

(H^-ihdt) / eT^^ {bo + hbi)deipn{7])dri 
The equation for the second order symbol: 



h^ 

2m ,/ron+fc 



eT^^ A^bideiph{r])dT] + 0{h° 



dtb2 + V^.S V^62 + r^A^5 62 = i^A^bi, 
2m 2m 



[ 62(0, X, r/, 0) =0. 



implies 

{H^ - ihdt) 



eft' 



'(60 + hbi + h%)de0n{v)dv = f ei^ A^b2de^n{v)dyi + 0(/i~) 

2m J^n+k 
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Therefore we can deal with functions bj, j > 1 fulfihing the recurrent equations 



1 i 

dfbj + V^5 V^6j + —A^S bj = —A^bj-i, 
2m 2m 



(4.5) 



, bj{0,x,r],e) = 0. 

Each solution bj is a symbol as in Definition 13.21 and therefore (thanks to Th. 13. 5p it defines 
a bounded operator Bj{t). The same holds true for the remainder operators: 



nj{t)ip = {2-Khy 



2ri + fe 



where r,- := A^bj. In order to prove it, we need the following 

2m 

Lemma 4.1. For all j > 1 the solution of equation ( [^.Jp fulfills the estimates: 
\bj\, \A^bj\{t,x,r],e) < I 

Cr{T)X-^{x,v)e-\'\\ {x,v,e)eTs. 
Proof We consider the following problem 

^C{t, X, 7], 9){t) = V^S{t, Cit, x, rj, e){T),7], 9) 



(4.6) 



(4.7) 



and define 



<^{T,x,r],9) :=exp|-^^ A^Sit,C{t,x,r,,9)ir),rj,9)dr^ 

in order to apply the theory of characteristics {{9,rj) fixed) and find the solution: 

i /■* 

bj{t,x,ri,9) = — <^{t-T,x,ri,9)A^bj-i{T,C{t,x,7],9){T),v,9)dT (4.8) 

By the iteration of this map, we have a direct linear relationship between Axbo and bj. Now 
we recall the estimates on 5o proved in Theorem 12.221 

\d^bo{t,x,r],9)\ < C+(r)e^"(^)^(^'^)e-l^l' y{x,r],9) e R^''+'' 

and the explicit analytic structure of S studied in Theorem 12.131 



S 



+ 9{t,x,9) 



(4.9) 



which implies the exponential behaviour of ^. The exponential upper bound for \bj\ and lAjj^jl 
follows directly from that. In the region Tg we recall the upper bound of type A~"(x, ry)e~l^l^ 



37 



we required for Axbo and the expansion AxS{z) = AxS{z) + G{z) with G £ (see Th. 
I2.16p . By using (|4.8p we obtain this second estimate also for \bj\ and lA^-ftjl. □ 

As a consequence, we can apply the boundedness result of Theorem 13.51 and state the 
existence of constants Kj{T) > such that ||7^j(t)|| < Kj{T)h'^^^ . By well known arguments 
related to the Duhamel formula we obtain the estimate: 

U{t) -Y,B,{t) <T ll^iv(s)|| ds < TKN{T)h^^\ t G [0,r]. 



j=0 



□ 



Now we clarify the relationship between the construction of the previous theorem and 
Chazarain's formulation |Ch) . as well as with the integral representation of Fujiwara [Fuj . 



Theorem 4.2. Let t S [0, to]; with tQ so small that the solution of the Hamilton-, J acobi 
equation does not develop caustics. Consider the construction of Theorem 1.1, truncated at 
any finite order J : 



J J 



VSj(t)(^:= V(27r;i)-" / / / e*(^(*'^'^'^)-<^'^>);i^6j(t,x,??,e) dr/ ^(y) (4.10) 



3=0 j=0 
Then: 

1. J J 

^ B,{t)^ = ^ Utit)^ + 0{h'+') (4.11) 

jr=0 «=0 

Here U^{t) is the term of order h" of Chazarain's FIO (ICh!^ ). 

2. 

J J 

B,{t)^ = U^it)^ + 0{h'+') (4.12) 

jr=0 0=0 

where this time U^{t) is the term of order of Fujiwara' s integral operator (JFu!]). 

Proof. In order to prove the first assertion, the main idea is to apply the stationary phase 
theorem to the oscillatory integrals ()4.10p with respect to 0- variables. In the same way, if 
we consider the stationarity argument with respect to (0, ry)-variables we obtain the second 
assertion. 

In the small time regime t G [0, to] there exists a unique smooth and global critical point 
6*{t,x,r]), solution of = V0S{t,x,r],6). This fact suggests us to consider the translated 
phase function around this point S(t,x,ri,9 + 6*(t,x,r])) with 9 £ -Bi(O) and symbol bo (see 
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Theorem 1.1) for which we choose the regularizing part as p{6) := {volBi{0))^^Xi{9), a C°° 
cut off function for the ball i?i(0). The compact behaviour of bj on the ^-variables follows as 
a consequence. The uniqueness of 9* and the compact setting in the oscillatory integral allow 
us to apply the stationary phase theorem to each integral in the ^-variables 

Bj{t,x,r])= [ eTi^^^'''''''^h^bj{t,x,r],9) d9 

obtaining 

Bj{t,x,T^) = e:^^(*'^''''^*VetV^'5(t,x,r/,r(t,x,r?))r^ e^''Wbj{t,x,i],e'') + 0{W+^) (4.13) 

where a = sgnVgS{t, x, tj, 9*{t, x, r])) and we have omitted (to simplify the exposition) the ex- 
plicit form of the higher orders symbols. Now we remark that the function S{t, x, rj, 9*) equals 
the phase used in the Chazarain's paper (the action functional evaluated on the classical curve 
with boundary conditions x and rj). Hence, by the uniqueness of the symbol expansion of the 
propagator in powers of ^ , we get the corrispondence between the symbols obtained as in 
(|4.13p and the ones obtained in the above-mentioned paper. This implies the equivalence of 
the two series (j4.1ip up to an order o{h'^~^^). By the same argument, applied this time to the 
integrals over u := {9,7]) and ^{t, x,y,u) := S{t,x,r],9) — {y,r]) 

Bj{t,x,y)= f [ eT^'^^''^''^'^'>-^y'''^h^bj{t,x,rj,9) d9 di]= [ e**(*'^'^'");i^6j(t, x, u)dn 

we use the uniqueness of the critical point u*{t,x,y). to get 

Bj{t, x, y) = eS*(*'^'2^'"*)n^'|detVi«>(t, x, y, u*{t, x, y))\'h'f'' bj{t, x, u^t, x, y)) + 0{h^+^) 

The phase function ^{t,x,y,u*) is the same used by Fujiwara and therefore also (|4.12p is 
proved. This concludes the proof of the Theorem. □ 



5 Multivalued WKB semiclassical approximation 

In this final section we prove Theorem 1.2, mainly applying the Stationary Phase theorem to 
the global FIO (|4.2p . in order to get a multivalued WKB semiclassical approximation of the 
Schrodinger evolution operator. 
Proof of Theorem 1.2 

We start by recalling that (as proved in Theorem 1.1) the ^-Fourier Integral Operator 

Bo{t)^:={27rh)-'^ f f f e-^^^'^''''''^'^^-^y'''^^bo{t,x,r,,9) d9 dr] ^{y) dy (5.1) 

JR" JR" JRfe 
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is a semiclassical approximation of the Schrodinger propagator for all t G [0,?"]. Under the 
particular hypothesis 

V{x) = l-\x\^ + Vo{x), sup ||VVo(x)|| < 1, (2t + 1)^ (r GN), 

we proved (see Theorems l2.18ir2.19l and l2.20p that the phase function has isolated and finitely 
many critical points; precisely the equation 

VeS{t,x,r],9)=0 (5.2) 

is solved on a finite open partition (x, r/) G M^" = |^ in such a way that on each Di there 

e=i 

are exactly i smooth solutions 6*{t,x,r]), 1 < a < i. This property allows us to apply the 
Stationary Phase Theorem (see [Ho2j vol. I) to the oscillatory integral in ()5.ip . The result is: 

Bo{t,x,r,) = [ ei^^'^^^'^'%o{t,x,ii,e) dO 

= ^et^(*'^'^'^S(*'^'^))|detVi5(t,x,r?,0^(f,x,r?))|-leX-'^6o(t,x,r?,C(i,x,^ 

a=l 

+ o{h) 

where cJq, = sgnVgS{t, x,r],9^{t, x,r])). In the small time regime t G [0, to] and for potentials 
V quadratic at infinity, it is well known (see i.e. (We]) that the graph of the Hamiltonian flow 

At := { (y , r/; X, p) G T^M" x T^M" \ {x,p) = 4>i^{y,r,)} 

= {iy,v;x,p) G ™" X ™" I p = V,5, y = VriS, = VeS} 

is globally transverse to the base manifold [x, rj) G M^", so the equation ()5.2p admits a unique 
global smooth solution 0*{t,x,ri). This simplified setting yields: 

Bo{t, X, 1]) = e|5(M'fl,e*(t,xfl)) |detV^5(t, x, t], e\t, x, r]))\-h e^- bo{t, x, r], 9*{t, x, rj)) 

+ o{h) 

which is the usual WKB construction, local in time. □ 
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